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I. INTRODUCTION 

The amplitude analysis of non-leptonic three body B 
', decays became an important tool to determine the CKM 
phases [H, S| , and also to observe CP-asymmetry. Us- 
ing this method, the Belle and BaBar collaborations have 
recently [J, [H] extracted a fraction asymmetry for the 
channel K^p'^. With this kind of analysis one 

could also explore the asymmetry associated with the 
interference between two neighbour resonances decaying 
into the same three body final state. 

The Dalitz plot analysis needs some a priori model, 
I with all possible dynamical components and a correct 
functional form to be used to fit the Dalitz plot distribu- 
tion. Arbitrary distribution functions can be used to get 
a good fit, but they have no physical meaning. 

The non-resonant component which, in general, is 
spread all over the phase space, can mimic other dynam- 
ical components, through the interferences with the reso- 
nances present in the same phase space. This shadowing 
phenomenon was observed in charm three body decays 
in the E791 experiment [1, 0], where the overestimated 
contribution of the non-resonant amplitude replaced, in 
a wrong way, the contribution (and the very existence) 
of the scalar mesons a and k. 

Belle has proposed, for the amplitude analysis of the 
' process — > K^Tr~^Tr~ a parametrization for the 
non-rcsonant amplitude given by: 

; Anr{K^7r^TT^) = a';:''e-°"''e'^"" +aYe-°"'''>e'^^\ (1) 

where oi, Si, 02 and 62 are the fit parameters and 
si3 = M'^{K^TT^) and S23 = AP{TT+n~) are the Dalitz 
variables. Originally the parametrization was also func- 
tion of the variable S12 — M'^{K^n^), but the term con- 
taining this variable turned out to give an insignificant 
contribution. 

The function above could fit the data with an accept- 
able confidence level around 1%. Certainly this distribu- 
tion was more useful for the Belle analysis than the usual 
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constant non-resonant distribution. However, it has to be 
employed with caution. 

The use of an empirical parametrization, with a form 
without dynamical content, may hide the physical mean- 
ing of the Dalitz plot and, even worse, may yield inad- 
equate parameters for the contribution of the resonance 
amplitudes and for the CP asymmetries. This is so be- 
cause in the Dalitz distribution the parameters are highly 
correlated among themselves and also with the non- 
resonant amplitude. The parametrization ([1]) was pro- 
posed to describe non-uniformities in the non-resonant 
three-body decays. Data suggest the formation of a tt~ 
with momentum predominantly smaller than the momen- 
tum of the and 7r+ in the process B^ /T+tt+tt^. 
How can we understand this? 

In principle there are two different cathegories of dy- 
namical processes contributing to the non-resonant struc- 
tures in the Dalitz plot. One is a genuine three-body 
decay that is a consequence of the partonic structure of 
the weak decay [1, |^ and therefore will be called par- 
ton interaction, or PI. The other comes from final state 
hadron-hadron interactions [l^, [O, [H, [3, and 
we shall call it FSI. In both approaches one is forced to 
make non-trivial calculation assumptions and the param- 
eters and form factors are not well known. However as we 
are dealing with Dalitz distributions, the comparison of 
these different contributions with the experimental distri- 
butions might discriminate which dynamical mechanism 
is more realistic to describe data. 

In Fig. 1 we show the diagrams relevant for the par- 
tonic description of this process. Figs, la and lb show 
the direct and penguin contribution to the "current in- 
duced" processes. Figs. Ic), Id) and le) show the direct 
(Ic)) and and penguin (Id) and le)) contributions to the 
"transition" processes and and Figs. If) and Ig) show 
the direct and penguin "annihilation" processes. 

Looking at diagrams la) to le) we can see that in the 
beginning both b and u quarks can carry a large momen- 
tum. When they emit a, W or a gluon, the bosons can 
also carry a large momentum and therefore, with the ex- 
ception of the d and d quarks (which come always from 
the vacuum and are soft) any of the quarks in the final 
state may carry a large momentum. This large momen- 
tum is then transferred to the final mesons. The final 
momentum distribution of the three mesons will be even- 
tually non-uniform but there is no reason for producing 
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Parton description of the weak decay 



a softer tt~ . The "entanglement" of the b, u and bosonic 
hnes will distribute more or less democratically the ini- 
tial high momenta. In sharp contrast to this situation 
are diagrams If) and Ig), where after the boson emission 
and absorption the high momenta are only with the u 
and s quarks. In the middle of the diagrams we see the 
7r~ made of u and d quarks pairs taken from the vacuum. 
This 7r~ will be comparatively softer than the negative 
pions of diagrams la) - le). 

When the final three mesons are produced through fi- 
nal state interactions we can also understand the softer 
7r~ in a very simple way: first the i?+ decays into 
and TT*' and then these interact creating if^, 7r+ and vr^. 
The first two mesons carry the valence quarks of the more 
energetic intermediate if + and 7r° and therefore will have 
higher momenta. The 7r~ is produced with quarks taken 
from the vacuum and therefore is softer. This is illus- 
trated in FigO Thus the appearance of a comparatively 
soft Ti~ is nothing but the manifestation of the "leading 
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FIG. 2: Born diagram for K+i^^ 
hadronic interaction. 
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particle effect" , so frequently observed in other processes 
in hadron physics. 

To the best of our knowledge, this kind of process has 
been first observed in a Dalitz plot by the Aachen-Berlin- 
CERN Collaboration in the reaction 7r~p Tr+Tr^n [l6l | 
and also in 7r~p K^K^n Shortly after the ex- 

perimental observation it was theoretically understood 
in [Tsj . In that paper it was explained in terms of double 
Regge graphs. This mechanism implies the production of 
a soft pion in the central region, responsible for a bump 
at the corner of the Dalitz plot, hence the name cornering 
effect. 

The analysis performed in [l^, [lO], suggests that in 
the chiral limit approximation, i.e. when m\ >> 
diagram If) vanishes. It also suggests that the most im- 
portant process is the one shown in Fig. Ig), which has a 
striking similarity with the diagram of Fig. [51 Both have 
the same final quark flow and final meson formation. 

From the above qualitative discussion and from figures 
Ig) and[2]we conclude that in the decay X+tt+tt" 
we expect to see a softer ■k~ which will yield a non- 
uniform Dalitz distribution. This effect comes both from 
the partonic weak decay and from final state hadronic in- 
teractions and has, in both cases, the same physical ori- 
gin: valence quarks form hard mesons in the final state 
while softer mesons are produced from the vacuum. Ne- 
glecting the flavor change (h — > s) in Fig. Ig), we can 
say that in both cases, PI and FSI, we have a leading 
particle effect. 

In what follows we will implement these ideas in a more 
quantitative way. In view of the uncertainties in the 
partonic description, specially the use of the factoriza- 
tion hypothesis, we will develop here only the hadronic 
final state interaction approach, for which the relevant 
lagrangian densities and form factors are better known 
from a large body of phenomenological studies at low 
and intermediate energies. Therefore we shall suppose 
that the non-resonant three-body decay K^tt~^tt~ 
proceeds through a hadronic scattering between two in- 



3 



termediate mesons and tt" in the two-body decay 



II. HADRONIC FINAL STATE INTERACTIONS 

Hadronic final state interactions have a long story and 
follow two different approaches: Regge models and me- 
son exchange models. The former have been considered 
in many works jTol . [l^ whereas the latter were discussed 
in [H, 13, [IB] ■ Regge theory is formulated in the high 
energy limit and in it all the amplitudes respect unitarity 
constraints. The interactions are represented by reggeon 
exchanges, the Pomeron being usually the most impor- 
tant one. As one moves to the low energy region, sec- 
ondary reggeons become important, additional assump- 
tions have to be made and the theory looses predictive 
power. It is not clear when the energy starts to be "high" 
and Regge theory has been applied to energies of the 
order of a few GeV. Meson exchange models describe 
hadronic interactions fairly well at intermediate energies 
(a few hundreds MeV). They contain uncertainties as- 
sociated with higher order diagrams, multiple exchange 
terms, coupling constants and the spatial extension of 
hadrons. These uncertainties are translated into form 
factors, which, in turn, can be either calculated (e.g. with 
QCD sum rules [2ll |) or simply parametrized and fitted 
to data. A positive aspect of meson exchange models is 
that one can use effective Lagrangians for them and so 
enforce chiral symmetry, which is known to be very im- 
portant at intermediate energies. On the other hand, in 
this kind of model, the amplitudes are not unitary. The 
lack of unitarity is believed to become a serious problem 
at increasing energies. These considerations suggest that 
with inceasing energies we should change the dynamical 
description, going from meson exchange to reggeon ex- 
change, but is not clear at which energy this should be 
done. In this work we shall study final state interactions 
with meson dynamics. 

In principle the B meson can decay into inumerous 
hadronic two-body intermediate states which will sub- 
sequently decay into X+tt+tt". These states include: 
K°TT+, K+n^, ri'K+, 7fK*+, r]K+ , r]K*+ , ujK+ , ujK*+, 
a+K°, a°K+, K*°tt+ , K*+n°... These mesons can un- 
dergo a two-to-three reaction, exchanging virtual mesons 
and creating the 7t~ . This is illustrated in Fig. [5] for the 
K~^TT^ intermediate state, which we will use as working 
example. A more careful analysis reveals that some of 
the above mentioned two-body states can be excluded 
because their interactions yielding a K'^n~^'!T~ final state 
would violate G parity (this is the case of the A'°7r+ and 
QqK'^). Some other states are much less abundant as 
two body decays. This is the case of the ujK'^ and ujK*^ 
states, which branching fraction is one order of magni- 
tude smaller than the K'^tt'^ one. Other diagrams are 
suppressed by dual topological unitarity (DTU) [l^l . Fi- 
nally, some other states, inspite of being different, are 
expected to yield amplitudes which are qualitatively very 



similar. This is, for example, the case of the states /C+tt^ 
and r]'K~^ . Replacing 7r° by rj' in Fig. [5]will lead to a sim- 
ilar diagram with the vertex rj' pn instead of tt pn. This 
will change the corresponding coupling constant and the 
form factor. The former is unable to distort the Dalitz 
plot and irrelevant for our discussion. In both vertices the 
form factor will suppress highly virtual internal p mesons. 
We therefore expect these two diagrams to give similar 
results. Even after a careful scrutiny and the elimination 
of suppressed diagrams and of "double counting" , there 
will remain a few diagrams to be considered. Here we 
shall work out in detail only the process depicted in Fig. 
[2l This will show schematically how the final state meson 
dynamics will lead to inhomogeneities in the Dalitz plot. 

Before starting the evaluation of the Feynman diagram 
given in Fig. [21 one last remark is in order. In other FSI 
calculations, such as [l4| the starting point is the B me- 
son, which decays into two off-shell mesons, which, in 
turn exchange a third virtual meson. This process is de- 
scribed by a loop diagram, from which the absorptive 
part is considered. In our Fig. [2]this procedure would be 
equivalent to close the lower part of the diagram forming 
a loop with a virtual K'^ and a virtual tt*^. Here, for sim- 
plicity, we take them to be real. In this way we will give 
emphasis to the creation of the soft 7r~ and consequent 
Dalitz plot distortion. The loop integral will introduce 
some smearing in our result and presumably change the 
normalization. Since our purpose in this work is to dis- 
cuss this process qualitatively, we postpone to the future 
a more general calculation, including effects of the kaon 
and pion off-shellness. 

The effective Lagrangians relevant to the calculation 
can be constructed from a chiral Lagrangian [2^. They 
are: 

^piT-K = gp-nivP^ ■ (tT X 9^7?), 

-Cpc^^ = gpu-Ke°''^^'^di3UJaPx ■ d^n. (2) 

With the above interaction Lagrangians, we can write 
the amplitude for the diagram in Fig. [5] as: 



A = 



(3) 



{{Pi - Ps)^ - ml) {{p5 - pa)^ - "^2) 

where pi and p2 denote the momenta of the and tt*' in 
the initial state; p3, p^ and p^ are those of the , tt^ 
and 7r+ in the final state of the diagram in Fig. (2] It 
is important to notice that pi + P2 = Pb, where ps is 
the momentum of the B meson, since we are considering 
the final state interactions in the two-body decay 5+ — > 

In order to take into account the effects of hadron in- 
ternal structure we follow ref. [l^l and introduce, in the 
amplitude, the form factor: 

. / A* \ / A-* 



+ [qi - mlY 



(4) 
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where and qp are the four momenta of the intermediate 
off-shell vector mesons, i.e., q^ = pi —p^ and Qp = p^ —p2 
and A is a cut-off parameter taken to be A = f GeV. 

Parametrizing the amplitude of the weak decay B"*" — *■ 
K~^TT^ through the weak decay coupling Gbktt, we can 
write the amplitude for the three-body decay, — )■ 
if+TT+TT^, in terms of the K~^tt~ {pk-k ~ Ps + Pi) and 
TT+TT^ (p^TT = P4 + Pb) momenta as: 



A 



-32iA^G 



K-KTT 



BKn9ujKK9p7r7T9puj7r 



((?2-4m2)(A4 + (g2/4-rn2)2) 



{P-K-K)aPBI3P\{PK-K)c 



{ql~Aml) (A4 + (g2/4-m2)2)' 



(5) 
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(15) 



(16) 



(17) 



(18) 



where qi ^ P ~ ps + 2p^Tr, q2 = P + Pb - 2pK-K and 
P = Pi — P2- In terms of these momenta, the three-body 
decay rate is given by: 



P-K-K -PK-JT 



(19) 



dV 



l-^i^TTTrl 
247r5TOB 



5 {{pnTv +PK-K -PbT - ml) S {{pi 



- PK-nT - m-l) S {{pB - P-mx f - m\) d'^p^.T, d'^pK-n- (6) 

Evaluating Eq. ^ and using the delta functions to 
perform some of the integrals, we finally write the decay 
rate in the standard form for the Dalitz plot: 



dsi3ds23 
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where, as in Eq. (P), S13 = pj^^, and S23 = pItt- In the 
rest frame of the B meson and due to the delta functions, 
in evaluating |-4i<-7r;T| one has to use: 



P^ = 2ml 



2m^ — m? 
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FIG. 3: The Dalitz plot for the non-resonant -> K+tt+tt" 
decay. 

For a given value of S13 ((m^r -I- mn)^ < S13 < (ms — 
TO^)^), the range of S23 is determined by imposing — 1 < 
cos ^ < 1 and is given by [25] : 



where X{x, y, z) = x'^ + y'^ + z'^ — 2xy — 2xz — 2yz. For 
the other momenta and scalar products we have: 
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Ef - ml - JEf - ml 



ql = 2si3-2mg + 2ml + Amj^-4P.pKTr, (H) where 



P.PKTT = EpEktt - \P\\pk-k\ cosOktt, 
P.pn-^ = EpE^r-^ ~ \P\\pt,^\{cos6k^cos9 



cos — 



{mji - ml -m% + S13 + S23 + 
{si3 - ml){s23 - mj.) 



mr 



(12) 



(13) 



(14) 



E* 



E* = 



si3 - + ml 



2pK-K 
i| - Si3 - 



2pK-n 



Before presenting our numerical results we would like 
to emphasize that we do not wish to reproduce the abso- 
lute normalization of the decay rate. The only purpose 
of our calculation is to show that the meson exchange 
mechanism is able to produce distortions in the Dalitz 
plot of the B+ K~^Tr~^Tr~ decay. 
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In Fig. [3] we show the Dahtz plot of a Monte Carlo sim- 
ulation of Eq. ([7]) . The distribution shown in the figure is 
consistent with the parametrization in Eq. (jlj but it goes 
to zero near the threshold. This result gives support to 
the expectation presented in section II: the K'^ and 7r+ 
carry the valence quarks from the intermediate K'^ and 
TT^ and are therefore hard, while the is created from 
the vacuum and is therefore softer. 

The Dalitz distribution of Fig. [3] can be parametrized 

as 



\AnriK~^Tr 7r+)p cx V^23^/i(s23) /2(si3) e' 



-Dsi 



where 



1 



1 + e i'-c-Pi)] 



(20) 



(21) 



Pi 



with D = 1.3232 x lO^^ GeV-^ a = 0.65 GeV-^ 
18 GeV2, C2 = 0.55 GeV-^ and p2 = 15 GeV^. 

The meson exchange mechanism considered above 
could play the same role in the Z?+ K^tt^tt^ three- 
body decay, where one does not see any inhomogeneity 
in the non-resonant amplitude. In the case of the , 
the three-body decay could proceed through final state 
interactions between the Ks and 7r+ mesons in the two- 
body decay Kgir^. In Fig. 2] we show the Feyn- 
man diagram for the —^ K~'k~^tt^ decay through 
the Kg 7r+ interaction. We see that there are no mesons 
that could be exchanged satisfying all the conservation 
laws required by strong interactions. This means that 
the diagram does not exist and the three-body decay 
Z)+ K~TT^-K^ does not proceed through final state 
interactions. A similar situation occurs in the thre-body 
decay i?+ I?^7r+7r+, since there is no possiWe me- 
son exchange that would lead the two mesons I?" and 
Tr"*" in the two-body decay D^n^ into the final 
three mesons: _D^7r+7r+. On the other hand, the three- 
body decay K^K~K^ can proceed through the 
two-body decay K~^7r~, followed by final state in- 
teractions between and tt^, as can be seen by the 
Feynman diagrams in Fig. [5l 

Therefore, in the case of the three-body decay -6° 
K'^ K~ , one expects some inhomogeneity in the non- 
resonant amplitude, consistent with a hard ^ a hard 
and a softer K~ . This effect was not included in the 
data analysis performed by the BaBar Colaboration in 
ref. (29j . where a homogeneous parametrization for the 
non-resonant component was used. 

To close this section we would like to comment the 
Belle results on baryonic three body decays. In this con- 
text the FSI approach was first used in [30| to study rare 
baryonic B decays. 

The study of the decays B^ ppK^ , B^ ppir^ 
(3l| revealed a low mass enhancement of the system p—p. 
However this enhancement in the low mass (close to 
the threshold) baryon-antibaryon system was not seen 
in the decays B ppj/ip and B" ApJ/ip [33 |. 
In terms of FSI we can very easily understand these 
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FIG. 4: Born diagram for KsTt'^ 
hadronic interaction. 
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K 7r+7r+ final state 



FIG. 5: Born diagram for K^tv' 
hadronic interaction. 
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results. The mass of the baryon-antibaryon system is 
low because the p is soft. In the case of light meson 
production we have the following sequence of decays: 
i3+ 7r° K'^pp. The intermediate pion splits 

into a proton and an antiproton and the latter interacts 
with the kaon exchanging a p or an w. The emerging p 
has all the antiquarks coming from the vacuum and will 
therefore carry low momentum. In the case of J/tp pro- 
duction the B goes first to an intermediate state with 
J/ip, such as, for example, J/ipT^^- The pion could then 
split into p and p but now the subsequent elastic interac- 
tion between the J/tp and the p is OZI suppressed, since 
the two mesons have no quarks in common. Therefore, 
in this case there is no soft p from FSI and no low mass 
enhancement in the pp system. This same problem has 
been addressed in |l(| where a similar conclusion was 
obtained. 
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III. CONCLUDING REMARKS 

We hope to have made clear that, for experimental 
purposes, it is highly desirable to have a parametriza- 
tion for the amplitude of non-resonant three body de- 
cays. This parametrization should have a physical basis 
and should be derived from theory. The already exist- 
ing parametrization given in Eq. ([1]) is not satisfactory 
from this point of view. In order to improve it we have 
considered one possible mechanism responsible for inho- 
mogcncitics in the Dalitz plot, which we have called final 
state interactions (FSI). We have sketched calculations 



of the Dalitz distributions with this mechanism. These 
calculations can certainly be improved, but already at 
this stage they can teach us how to obtain distortions in 
the non-resonant three body decays. FSI trace back the 
observed Dalitz distributions to a manifestation of the 
leading particle effect. 
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